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ABSTRACT 


The  use  and  processing  of  the  FM/CW  signal  for  radar  and  acoustic 
sounder  systems  are  examined  In  this  note.  This  signal— along  with  real¬ 
time  digital  processing  via  minicomputers— Is  currently  being  used  by 
several  groups  for  HF  over-the-horizon  radars.  A  comparative  analysis 
of  the  different  processing  techniques  for  general  radar  applications 
has  yet  to  be  undertaken.  This  note  therefore  attempts  to  promulgate 
details  of  these  techniques  so  that  they  may  find  use  In  other  systems. 

An  example  Involving  an  HF  backscatter  radar  Is  used  to  permit  the  reader 
to  see  how  the  techniques  are  applied  to  an  actual  problem. 

A  linearly  swept-frequency  signal  format  Is  used  In  a  100%  duty- 
factor  mode.  In  the  receiver,  a  replica  of  the  linear  FM  signal  Is  mixed 
with  the  received  waveform  at  an  offset  such  that  the  desired  range  window 
is  observed  with  the  lowest  possible  IF  frequency  variation.  This  pulse 
train  Is  then  analog-to-digital  (A/D)  converted  and  ready  for  computer 
processing.  Two  techniques  are  described  and  analyzed  for  digitally 
processing  the  signal  via  the  Fast-Fourler-Transform  (FFT)  algorithm. 

The  first  Is  a  double-FFT  process;  the  first  FFT  set  Is  done  within  a 
pulse-repetition-interval  (PRI)  to  give  range  Information.  The  next  FFT 
set  Is  done  over  N  PRIs  to  give  Doppler  Information.  In  the  second 
technique,  a  single  long  FFT  Is  used  over  N  PRIs,  giving  simultaneously 
both  range  and  Doppler  information.  It  Is  shown  that  both  techniques  are 
Identical,  In  that  they  produce  the  same  Information  and  require  the  same 
number  of  computer  steps  In  executing  the  required  FFTs.  Both  techniques 
yield  unambiguous  range  and  Doppler,  for  both  discrete  and  distributed 
targets;  the  note  shows  how  and  where  this  Information  Is  contained  In 
the  processor  output.  The  note  also  describes  how  two  weighting  functions 
are  normally  applied  to  the  pulse  train  time  samples  to  reduce  objection¬ 
able  range  and  Doppler  sidelobes.  Finally,  simple  "cookbook"  rules  are 
given  for  obtaining  the  signal  and  processing  parameters  based  on  the  radar 
and  target  range/ velocity  specifications. 
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FM/CW  RADAR  SIGNALS  AND  DIGITAL  PROCESSING 
Donald  E.  Barrick 

1.  OBJECTIVE 

The  objective  of  this  note  Is  to  present  a  simple  and  concise 
analysis— backed  by  an  example— of  the  application  of  an  FM/CW  signal 
format  In  radar  systems.  It  Is  shown  how  both  time-delay  (range)  and 
Doppler  (radial  velocity)  Information  can  be  extracted  unambiguously. 

2.  APPLICATION 

For  the  sake  of  Illustration  throughout  these  notes,  we  pick  the 
following  application  and  example.  The  HF  radar  carrier  frequency  Is  to 
be  10  MHz.  Sea  scatter  Is  to  be  observed  from  the  radar  out  to  a  range 
of  150  km  (corresponding  to  time  delays  up  to  1  millisecond  In  a  back- 
scatter  radar).  It  Is  known  that  HF  sea  scatter  Is  confined  spectrally 
to  frequencies  within  about  1/3  Hz  of  the  carrier.  Therefore  a  pulse- 
repetltl on-frequency,  f^,  of  1  per  second  Is  selected  so  that  all  echo 
Dopplers  within  +0.5  Hz  of  the  carrier  will  be  displayed  unambiguously. 
To  show  sufficient  detail,  a  Doppler  processing  resolution  better  than 
0.02  Hz  Is  desired,  and  a  range  resolution  of  the  order  of  1.5  km  Is 
desired;  the  latter  two  requirements  In  an  ordinary  pulse-Doppler  system 
translate  to  a  coherent  Integration  time  exceeding  50  seconds  and  a 
signal  bandwidth  of  100  kHz,  respectively, 

3.  TRANSMITTED  WAVEFORM 

We  select  a  100%  duty  factor  signal  whose  frequency  sweeps  upward, 
linearly,  over  one  pulse-repetitlon-interval  T^  (T^  =  1/f^  «  1  sec  for 
our  example).  Since  a  100  kHz  signal  bandwidth  Is  desired,  the  signal 
can  be  written 


Vy(t)  ■  cos[aj^t  +  n  Bf^t*]  H  cos[((»j(t)]  (1) 

for  -T^/2  <  t  <  T^/2.  It  Is  assumed  that  the  signal  Is  periodic,  and 
hence  phase-coherent  from  one  repetition  Interval  to  the  next. 


/ 


fr. W'.»TS  V*'V'  ’ 

r^ . — 


Since  the  instantaneous  frequency,  fy{t),  1s  the  derivative  of  the 
phase,  we  have 

1  d^rCt) 

^  /^\  _  I  I  _  ^  ,  r\^  A.  /a\ 


^T<‘)  ■  ST  -V  '  fc  ^ 


where  here  f^  ■  10  MHz,  f^  *  1  Hz,  and  B  =  100  kHz.  Thus  It  can  be 
seen  that  the  frequency  excursion  of  fj(t)  over  one  pulse-repetition 
Interval  Is 

Afj(t)  *»  B  *  100  kHz.  ( 


The  amplitude  of  the  transmitted  signal  Is  taken  to  be  unity.  The 
plot  of  signal  frequency  vs  time  Is  shown  In  figure  1. 

4.  RECEIVED  WAVEFORM 


Figure  1.  Frequency  ve  time  of  tranamitted  and  delayed/Doppler  shifted 
received  aignala. 
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The  received  signal  Is  both  delayed  In  time  and  shifted  In  Doppler. 

To  Illustrate  the  situation,  we  assume  that  we  have  a  discrete  target  at 
range  15  km  and  travelling  radially  away  from  the  radar  at  v>5  m/s  (e.g., 
an  ocean  wave).  At  time  t»0,  the  target  Is  exactly  at  Rj«15  km  from  the 
radar.  After  that,  Its  range  Is  a  function  of  time  as 

R(t)  -  R,+vt  .  (4) 

The  received  signal  from  this  discrete  target  Is  thus  just  a  replica 
of  the  transmitted  signal,  but  multiplied  In  amplitude  by  a  factor  A 
and  delayed  In  position  by  a  factor  t^,  where  t^j  »  2R(t)/c.  It  Is  thus 

Vp(t)  »  AVy(t-tj)  ■  Acos[a)j.(t-tjj)"hiBf^(t-t^j)*]  ;  (5) 

Its  frequency  Is  shown  In  figure  1  as  the  dashed  curve. 

5.  DECHIRPEO  SIGNAL 

Now  after  RF  amplification,  we  mix  the  received  signal  with  a  replica' 
of  the  transmitted  signal;  this  Is  represented  mathematically  by  sub¬ 
tracting  a  phase  <})j(t)  from  4>.j.(t-tj)  to  give  a  signal 

Vj(t)  ■  APy(t)cos[u)^(t-tj)  -  u)pt+TiBf^(t-tj)*  -  irBf^t*!  .  (6) 

There  Is  also  a  sum  signal  with  phase  4>i-(t)+((»y(t-tjj) ,  but  It  Is 
near  2w  (twice  the  carrier),  and  hence  removed  by  filtering.  The  function 

V 

Pj(t)  denotes  a  pulse  of  unity  amplitude  and  width  T,  where  here, 

T  » 

Thus  the  mixture  of  the  two  sawtooth  frequency  waveforms  and  their 
subtraction,  as  shown  In  figure  1,  produces  a  signal  whose  frequency 
format,  fj(t).  Is  as  shown  In  figure  2.  The  two  frequtncles  are 


and 
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The  Intermediate*  dechirped  signal  can  be  represented  as  the  sum 

of  two  pulse  trains  a.®  rhown  In  figure  2.  One,  v^(t),  Is  at  frequency 

f^,  and  the  width  of  these  pulses  Is  T  ■  T^-t^j.  The  other,  v^(t).  Is 

at  frequency  f^,  and  the  width  of  the  pulses  Is  T-t^j.  U  will  be 

possible  to  eliminate  v  (t)  by  filtering  If  f  »f  ;  such  will  be  the 

2  2  1 

case  here. 
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Figure  2.  Frequenoy  and  amplitude  plate  ve  time  of  reoeived  signal  after 
deohirping. 


Therefore,  we  are  left  with  a  single  pulse  train  to  analyze,  as 
represented  by  equation  (6).  It  Is  possible  to  re-certer  the  time  origin 
so  that  It  falls  In  the  middle  of  the  first  pulse;  this  Is  done  In 
figure  2.  The  frequency  and  phase  from  pulse  to  pulse  are  changing  very 
slightly,  however;  we  will  analyze  this  now. 

5.1  -T^/2  <  t  <  T  ./2 

Let  us  simplify  the  phase  in  the  first  pulse;  denote  Internal  t1me,t, 
within  this  pulse  as  t^.  Using  t=t^  and  tj»2R/c»2Rj/c+2vt/cHtp+2vt^/c, 
(where  tgH2Rj/c  Is  tiie  initial  delay  of  the  signal),  we  have 

4’i(ti )*<^T^^i ^ 

+  2»t-2  J 

f- (1  -  |)t|  .  (8) 

Thus  we  have  three  co'itrlbutlons  to  the  phase:  a  constant,  a  linear 
term  In  time  t^ ,  and  a  quadratic  term  In  time,  t|.  For  the  parameters 
of  the  example,  however,  the  quadratic  phase  term  Is  always  small  within 
the  Interval  -T^/2  <  t^  <  T^/2;  e.g.,  at  t^=T^/2,  it  is  of  the  order  of 
0.005  radian.  Also,  It  can  b'j  shown  that  the  second  term  In  the  linear 
factor  Is  small  compared  to  the  first  term  and  is  also  much  less  than  one 
radian.  Of  course.  In  all  cases  under  consideration  here,  v/c  «  1,  i.e., 
target  velocity  Is  small  compared  to  propagation  velocity.  Therefore  we  have 

<>j(t^)=.4',-2ir[  |£-f^.+Bf^t„]t^  ;  (9) 

hence  within  the  first  pulse,  the  frequency  fj  Is 

•  T 

As  can  be  seen,  this  frequency  offset  (also  shown  In  the  preceding 
figure)  consists  of  two  terms:  the  first  due  to  the  target  velocity  and 
the  second  due  to  the  time  delay  (or  range)  to  the  target  (tj»2Rj/c“n  1 


millisecond  for  R,»15  km).  For  the  exatiiple  selected  here,  the  second 
term  (range  term)  Is  larger;  I.e.,  Bf^tj*10  Hz.  Thus 

It  Is  not  possible  to  separate  range  from  target  velocity  by  measuring 
frequency  f ^  within  a  single  pulse. 

5.2  (2n-l)T^2  <  t  <  (2n+l)T^/2 

Here  we  want  to  examine  the  phase  In  the  n-th  pulse,  assuming  that 
the  n»0  pulse  Is  the  one  centered  at  t=0.  Again,  we  describe  the  time 
within  the  n-th  pulse  (from  Its  own  center)  as  t^.  The  time  delay  to 
the  target,  t^,  however  Is  now  given  by 

t^j  »  2R/C  =  2Ro/c+2vt/c  =  to+2v(nT^+t^ )/c  ,  (11) 

where  we  describe  time  to  the  center  of  the  n-th  pulse  as  nT^.  We  can 
now  substitute  this  into  the  phase: 

=  V»r‘d>  -  M*)* 

■  -  -ct.  -  T  -  “c  T  "V  *  T 

-  ZirBy  V  (nVt,)l  . 

After  expansion  and  elimination  of  terms  which  are  small  compared 
to  others  and  also  small  compared  to  one  radian,  we  have  (assume  n  £100) 

■  ♦.  ■  ^  ■  2”'  T  V®V.  *  T  '"’‘i  ‘ 

hence  the  frequency  in  the  n-th  pulse  Is  the  quantity  In  square  orackets, 

1  .e. , 


~  VBf^t  +  ^  Bn. 
c  c  roc 


(14) 


Comparison  of  (14)  with  (10)  shows  that  the  frequency  In  the  n-th 
pulse  Is  Identical  to  that  in  the  first  pulse,  with  the  exception  of  the 
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third  term.  The  explenetlon  for  the  third  term  \%  limpU.  It  merely 
means  that  the  target  Is  moving  from  pulse  to  pu1se»  and  Its  range  at  the 
center  of  the  n-th  pulse  Is  R,+c(2v/c)nT^/2  “  Rj+vnT^,  as  we  would 
expect.  Since  we  want  to  Integrate  over  as  many  as  100  pulses,  the  third 
term  Is  not  negligible  as  n  Increases;  e.g.,  at  n<*100,  ^  «  ^Hr. 

Two  other  effects  occur  within  the  pulse;  Its  width,  being  T''Tj,-t^ 
changes  very  slightly  from  pulse  to  pulse.  Since  T^,»l  second, 
t^j  ■  tj+(2v/c)nT^,  we  have  for  n^l,  T  •  I-IO"**  s;  for  n«100  we  have 
T  ■  I-IO’**-  j  X  10"*  s.  Thus  the  change  In  pulse  width  Is  negligible. 

A  very  Important  second  effect,  however.  Is  the  change  in  phase  from  pulse 
to  pulse,  as  represented  by  the  second  term  In  (13).  This  phase  change 
shall  In  fact  prove  to  be  the  basis  for  the  Doppler  processing.  As  stated 
earlier,  all  of  this  assumes  that  the  transmitted  signal  Is  phase-coherent, 
I.e.*  ♦j(t+T^)-4ij(t)  »  non-varying  constant. 

6.  DOUBLE-FFT  DIGITAL  PROCESSING 

Here  we  want  to  demonstrate  how  a  double  Fourler-transformatlon 
process  can  be  used— often  In  real  time  because  of  the  discovery  of  the 
digital  fast-Fourler-transform  (FiT)  algor1thm--to  produce  a  time-delay 
(range)  and  Doppler  (velocity)  display  of  the  radar  target  data*.  The 
first  Fourier  transform  process  Is  performed  over  a  pulse  repetition 
period,  T^  (I.e.,  within  a  pulse)  to  obtain  target  range.  The  second 
Fourier  transform  Is  then  performed  over  several  pulses  of  these  data  to 
obtain  target  Doppler  or  velocity. 

First,  let  us  perform  a  Fourier  transform  on  a  single  pulse.  This 
is  shown  in  figure  3.  We  have  a  pulse  of  width  T  =  amplitude  A, 

and  frequency  f^  given  by  (14).  To  perform  this  Fourier  transform 
digitally,  one  must  sample  the  pulse  M  times  within  the  time  period 
T^.  The  number  M  depends  upon  the  maximum  value  f^  can  assume,  and 
M/T„  must  be  at  least  twice  this  value,  I.e.,  2f,„„^,  according  to  the 
Nyquist  theorem.  For  the  problem  considered  earlier  where  we  want  to 

*Thi8  technique  is  currently  being  used  by  the  Stanford  Research  Institute 
for  real-time  processing  of  HF  ionospheric  radar  signals  at  their  Wide 
Aperture  Reseeurch  Facility  (WARF);  (Sveeney,  et  al. ,  1971). 
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Tho  Fourlor  tronsforw  of  tho  pul so  U  thon 

T/2 


V,.(f) 


V,„(f)  -j  Mco.»,„(t,  ^‘-1 

-T/2 

AT  (s1n[2n(f-fin)T/2l  -U,+12!rfc- 

^  l2n(f+fi„)T/2]  “■  ^ 
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Thit  PouHnr  triniforiN  1i  thoim  In  figure  3,  Slnnt  m  itirted  with 
{l^tOO  MlnlMuM)i  wt  obtiln  UMiplei  In  tht  frtgutncy 
doi»1n  fro*  ^*"^11**%  poiltivt  ¥tlut»  of  frt- 

ouincy,  Thou  itwpiti  iro  eonglox  In  MnofoK  01  ovldonctd  by  tht  OKponontlil 
phtit  foctor  cunUlnIng  ond  Iwf^^nT^.  Thu«  wt  concoptuilly  htvt 
M/I  r«ngt  bint  (M/l*100  hort)*  pontlttlng  ui  to  rtollio  tho  1,1  kRi  nngt 
roiolutlon  ovtr  t  ISO  bit  windowi  «i  InltUlly  ttlpulotod.  Noto  thot  ooc'h 
N/l  fotolutlon  tloitont  oftor  tho  firtt  FFT  cm  bo  conildorod  •  rongo  bln 
to  long  01  tho  Oopplor  tont»  (tv/c)fQ»  It  titoll  coitportd  to  tho  rongo 
tonn,  Bf^tii  thit  It  truo  for  tho  oxoiplo  contidorod  horo.  Sinco  ooch 
puUo  It  opproximouly  1/T  wido  ot  tho  holf-powor  point  (T<*T^  •  1  toe 
horo)»  wo  tihould  bo  oblo  to  rotolvo  100  torgott  In  rongo  bocouto  tho  width 
of  ooch  FFT  pulio  In  thit  100  Hi  window  It  1  H<.  Honco  oftor  ono  FFT 
procoit  within  0  pulto  wo  hovo  rongo  Informotlont  but  no  Dopplor  Infontotlom 
wo  turn  now  to  oxtroctlon  of  Dopplor, 

Noto  thot  If  wo  itort  with  tho  firtt  pulto  ot  n*l  ond  do  thit  FFT 
procott  on  ooch  pulto,  wo  obtoln  0  Fourlor  trontfonn  n  timot,  whoro  wo 
ovtumo  nsN  (tome  uppor  voluo),  Sinco  tho  troquoncy,  f,_i  ond  photo, 

2xf(.  ^  nTf.,  shift!  slightly  from  pulto  to  pulto  duo  to  torgot  volocity 
(os  givon  In  (14)),  this  s1n»/n  pulto  In  tho  froquoncy  domoln  will  chongo 
vory  slightly  oftor  ooch  Fourlor  tront formation.  Sinco  our  digital  FFT 
U  capoblo  of  producing  numbors  at  H/2  discroto  points,  (IS)  should  roolly 
bo  writton  with  f  roplacod  by  <  "  ^  -M/2  sms  M/2. 

Thus  tho  first  FFT  process  on  M  samplos  within  a  pulso  gives  M/2 
range  bins  for  each  pulse.  For  each  successive  pulse,  this  FFT  gives  M/2 
additional  positive  froquoncy  samplos.  Digitally,  we  store  each  M/2 
samplos  In  rows  of  0  matrix,  as  shown  In  figuro  4,  until  wo  have  N  rows. 

Thus,  wo  hovo  on  M/2-by-N  matrix  whoso  columns  so  for  roprosant  rongo 
bins. 

Now,  wo  perform  another  FFT  over  each  column,  or  range  bln.  This 
will  require  N  points  altogether.  Each  matrix  olomont  Is  a  complex 
number  whoso  value  changes  In  a  column  because  tho  froquoncy,  f^^i  and 
tho  phase,  2nf^*  ^  nT^,  are  changing  from  swoop  to  swoop.  Since  each 
of  tho  N  vortical  olomonts  comes  from  a  different  pulse  T^  sec  opart. 
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NT^  ite  irt  ftquiftd  to  fill  this  lootflN.  AliOi  n  cm  bt  folitod  to  time  from 
tho  fint  puUt  by  uio  pf  tmY^i  or  n"t/T^  (iQilni  1<  n<  N).  Htneo  tich 
colymn  It  rmlly  «  function  of  tlmoi  md  tho  N  eoluitn  oloittnti  cm  Ot 
eonildtrod  (dlQlttl)  imploi  of  thli  tlim  function. 

To  Fourlor  tronifont  ovtr  •  typicol  column  (tty  tho  m-th)i  lot  ui 
ogoln  rofor  to  our  OMimplo  for  tho  Urgot  it  R|»16  kmi  this  tirgot 
will  oppoor  In  tho  m-lO  bln  for  K/S«100.  At  wo  tow  boforo*  this  producoi 
fji^  1*  10  ♦  ^  •  10’*  Hi,  Thui  for  n  running  from  1  to  lOO-- 

eorroipondlng  to  timi  running  botwoon  1  md  100  locondt—two  thingi  hippon 
to  tho  poll t1 VO  pulio  In  tho 


fift. 


m-th  rmgo  blnj  Iti  ompll- 

tudo  chongoi  illghtly  duo  to  j  .Kongo  “  “  P 

tho  ihlft  of  tho  i1n«/w  .  .  ^jm/A 

pulio  bocouit  of  f|g^t  ond  #  \ 

Iti  photo  chmgot.  Tho  I  S,*  ’  *  ‘  *  ^M/i  \ 

omplltudo  vorlotlon  from  n*!  I  ,  | 

to  n«100  It  slow.  For  tho  I  •  I 

OMomplo  g1von»  tho  thift  In  I  '  .  I 

tho  pul  to  duo  to  f,„  3*^*  I  m*  •  •  lun*  *  ‘  nM/i  I 

ovor  N«100  pul  tot  I  tho  3  dB  1  •  I 

width  of  tho  t1n«/«i  pulto  It  \  ’  / 

1/T.l  H<  ,*,11.  th.  toul  \  5  s,,.  .  .  S^.  .  .  S^/./ 

width  bttwMn  th.  flrtt  null.  '  “ 

1$  2/T"2  Hi.  Honco  tho  Figui^  4.  Matria  amtaining  renga-DoppUr 

,  ....  yntnib^rB  ohtaimd  with  doubU- 

omplltudo  vorlotlon  within  o  proocat, 

column  Is  slight*  ond  con  bo 

roprotontod  In  moot  costs  by  o  constont  or»  for  more  occurocy.  by  o  constont 
plus  0  smoll  lintorly  vorying  term;  tho  I'otults  will  not  diffor  slgnl- 
flcontly  for  olthor  cost.  Honco  wt  represent  tho  omplltudo  by  o  constant 
(1.0. .  s1nt2iT(f-f,^/2>^/2^/f2n(f-f^jj/2^^^2].  its  voluo  midway  down  the 
column  whoro  n«N/2)  ond  loovt  tho  second  rtprtsontotlon  os  on  exercise 
to  tho  Intorosted  rtador. 

Thus  tho  only  variation  now  within  tho  column  (at  tho  positive 
frequency  corresponding  to  m)  Is  the  phase  factor.  I.e.. 


$  .  ,  . 
II 

S  Its 

in 

^M/l 

$  ... 

It 

S  lit 

in 

^sM/i 

^nr  •  • 

S  %  %  % 
lun 

^nM/i 

^NI*  •  • 

S«  •  •  • 
(In 

^HM/i 

Figui^  4.  Matria  aemtaining  rungt-Dopplar 
ftumbtrt  obtained  with  doubU- 
FFT  proocat. 
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'mu 


.  K(,),'**^'  ¥  "V  .  K(f),'‘*'c  ¥  t„  . 


whtrt 


(18) 


In  th«  rlghtmoit  txprtsiloni  nT^  hai  btttn  r«pl«c«d  by  to  rtpmttnt  tho 
dUwhitt  flow  of  tlNW  from  pulio  to  pulit.  Tht  Fourltr  triniform  of  this 
quantity  ovtr  t^  from  0  to  NT^  Is 


^  fe)NV^l 


whort 


(17) 


htrt  iQiln  wt  should  noto  thst  our  dlqltnl  FFT  dots  not  rtslly  glvt  a 
continuous  variation  ovtr  f  (frtqutncy).  but  will  computo  valuta  at  N 
discrttt  frtqutncy  points.  Tht  qutstion  arista  at  to  how  wt  should  choost 
thtst  N  frtqutncy  points»  i.t.»  how  widt  a  frtqutncy  window  do  wt  want 
to  display.  Sinca  our  PRF»  f^(f^*l  Hz  htrt)  rtsults  in  an  unambiguous 
Doppltr  of  Hi,  wt  would  logically  stltct  J  f^(»  ^  Hi  htrt)  so 

as  to  display  all  of  tht  unambiguous  Doppltr  window.  Than  tht  frtqutncy 
window  in  Doppltr  will  bt  from  -fr,,.^  to  ’♦‘fc^ax  •  spacing  fn^^^/N , 
which  turns  out  to  tvtry  2f|^^/N  Hi,  or  1/100  Hi  htrt.  Nott  also  in 
(18)  that  if  Doppltr  shift,  txcttds  jr  fr*l/2T^»  thtn 

from  pulst  to  pulst  wt  will  bt  sampling  at  lass  than  the  required  Nyquist 
sampling  rata.  Hence  our  pulse-repetition  frequency  (PRF),  f,.,  must 
always  be  at  Ucmt  twice  as  great  as  the  maximum  expected  Doppler  frequency. 

Observe  now  an  important  fact  in  (17):  the  displacement  of  the 
sinx/«  pulse  resulting  from  the  second  Fourier  transformation  over  the 
columns  occurs  at  ^  f  •  This  is  precisely  the  Doppler  shift  that  results 
from  a  target  at  (radial)  velocity  v  with  a  backscatter  radar  having 
carrier  frequency  f^  Furthermore,  the  3  dB  width  of  the  pulse  represented 
by  (17)  is  1/NTp  Hi,  as  shown  in  figure  5.  Thus  we  produce  N  (or  100) 
Doppler  frequency  points  every  f^/N  Hz  (or  .01  Hz  here)  having  a  Doppler 
resolution  of  1/NTj.  Hz  (■  .0^  Hz  hero).  Since  NT^  is  the  coherent 
integration  time  (in  this  scheme.  It  is  the  time  required  to  fill  the 
matrix),  1/NT^  is  exactly  the  Doppler  resolution  one  would  expect  from 
any  coherent  pulse-Doppler  radar. 
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Thtrtfort,  In  \  havt  dont  two  ttti  of  PfTt.  Ont  tot 
within  oich  pulu  it  N  ^  nti  to  pivo  M/2  rtngi  blnti  thm  bint 
iro  tho  tloiMntt  of  •  row  of  •  MtHK.  Tho  itcond  itt  1i  ovtr  N,  pulton, 
or  ovor  tho  N  column  olomonti  of  tho  motrlM,  to  givo  N  Doppltr  bint 
for  oich  rongo  bln.  Noto  thit  tho  originti  tirgot  rtngo  olio  contolnod 
•  imoll  offiot  duo  to  Doppltr.  If  this  offiot  Is  objoctlonoblo,  It  con 
now  bo  rtmovtd— In  tho  cito  of  o  d1  scroti  Urgot»by  using  tho  Dopplo. 
Information  to  corroct  tho  torgot  rongo. 

A  llttlo  thought 
will  show  thot  this 
procott  olto  works 
for  dlstrlbutod  tor- 
gots  such  os  roln  or 
too  wovos.  If  ono 
hos  mony  torgots  In 
0  rongo  bln  (soy  L 
torgots),  ho  hos  I 
torms  In  (IS),  ond 
ooch  olomont  In  tho 

motrix  1$  roolly  ftgruro  S.  Doppltr  tptotrm  afttr  »«oond  tran§- 

tho  sum  of  L  such  formation  u>tthin  a  given  range  bin, 

torms.  Tho  second 

PFT  ovor  tho  columns,  thor*oforo  by  superposition,  gives  L  torms  In  (17); 
If  ooch  of  tho  L  scattorers  In  the  bln  (representing  the  distributed 
torgot  complex)  has  a  different  velocity,  then  the  resulting  Doppler 
spectrum  for  the  L  wurgets  will  consist  of  L-sIruc/x  pulses  at  different 
positions,  as  given  by  (17)  and  shown  In  figure  5.  Thus  a  continuous 
Doppler  spectrum  represented  by  the  sum  of  L  scatterers  with  many  differ¬ 
ent  velocities  and  scattering  amplitudes  will  result,  as  would  be  expected 
In  any  coherent  pulse-Ooppler  radar  system. 

7.  SINGLE-FFT  DIGITAL  PROCESSING 

Now  we  examine  another  technique  for  extracting  range  and  Doppler 
Information  from  the  same  signal.  This  Involves  a  single,  long  FFT  over 


.i  '“•»r?n.,»i.a^\^- wwvmu.s fi  ; .'r — t-tt^tv  tv.  •.^-  v-  ■  Tf-:mvn»T"— — 
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tht  tim  N  puitti.  TMi  xtchnlqut  1i  uttd  by  the  Rome  Air  Devtlopment 
Center  for  lome  of  Iti  Hf  over-the-hoHion  redert  (Eddy*  1973).  It 
Involvei  the  lenie  number  of  computer  operetloni  ei  thet  deicribed  In  the 
preceding  section. 

Here  we  will  drew  heavily  on  much  of  the  materiel  In  the  preceding 
section.  Since  we  have  a  maximum  frequency  f^^  In  our  pulse  T^(»T) 
seconds  loitg,  wc  rcqult's  M«2f,^*  Ty,  samples  per  pulsei  as  before. 
Performing  the  FF7  over  N  pulses  gives  a  total  of  M  x  N  samples  per 
transform.  Let  us  analytically  find  an  expression  for  the  Fourier  trans¬ 
form  of  this  pulse  train  first.  To  do  this*  we  can  use  superposition  to 
express  the  Fourier  transform  of  the  pulse  train  as  the  Fourier  transform 
of  each  pulse  as  though  It  were  all  alone: 

(an^DTy/a 


Vi(f) 


n-i  t 

E 

n«0  J 

(wn-l)Tj,/2 


'In 


(t).-'*”'*  dt 


(U) 


Here  we  reexpress  the  phase  appearing  In  as  given 

In  (13)— In  terms  of  continuous  time*  t»  rather  than  time  within  a  pulse* 
t^.  This  Is  done  by  substituting  t^*t-nT^  Into  (13)  to  give 

♦j^(t)  -  ♦,  +  2iT[Bto+  ^  BnT^]n-2TT[  ^  f,.  ♦Bf^t,+  ^  Bn)t  .  (19) 


Using  this  In  (18)  and  performing  the  Indicated  Integration,  we 
obtain  (with  the  approximation  T^xsT) 

,  .  -12nfnT« 

Vi(f)  -Y!  VT«(f)e  . 


-E  Vini 

nSt) 


(20) 


where  Vj^(f)  1$  given  In  (15)  anJ  discussed  In  that  section 

Now*  to  perform  the  sumnatlon*  we  make  the  same  assumptions  as 
before:  I.e.*  that  of  both  the  amplitude  and  phase  variations  over  n 
which  occur  In  Vj^(f),  only  the  phase  variation  Is  Important.  Also 
we  use  only  the  first  s1n«/a;  function  In  (15)  since  It  represents  positive 
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frtqutncittt  an  Identical  result  obtains  for  the  second  term  representing 
the  negative  frequencies,  Therefore.  (20)  becomes 

Vi(f)  -  K(f)V  ^  (21) 

ni^O 

where  we  have  used  (16).  The  above  summation  can  be  performed  by  using 
the  Identity; 


.1Na/2. 


sin 

$ 


(22) 


Thus  we  obtain 


s1n[2lKf-  ^  fc)NTr/2]  ^-12n(f-  ^  fc)(N-l )T^/2 
Sln(2it(f-  ^f^)T^/21 


(23) 


In  the  above  final  result,  the  complex  exponential  factor— as  well 
as  the  residual  phase  factor  e'*^°  contained  In  K(f)— Is  not  Important 
becausQ  1t  has  unity  amplitude.  However  we  note  that  the  s1nNx/s1n.T  Is 
much  like  the  s1nNx/N«  function.  It  contains  a  peak  at  x«0  and  side- 
lobes  away  from  the  peak;  It  Is,  however,  periodic  whereas  sInNx/Nx  1s 
not. 


Now,  we  note  that  the  FFT  does  not  actually  compute  a  continuous 
function,  Fj(f),  but  a  transform  at  MN  positive  and  negative  frequency 
points.  Since  the  maximum  frequency,  fmax*  Is  determined  primarily  by 
the  maximum  target  range  desired,  we  have  MN/2  poeitive  frequency  points, 
and  hence  a  value  of  Vj(f)  computed  every  f  =  2f ,^^y/(MN)  Hz  along  the 
positive  frequency  axis. 

To  see  how  a  discrete  target  will  appear,  we  plot  first  In  figure  6 
the  broad  function  representing  the  integration  over  a  single  pulse,  I.e., 


^  sin[2Tr(f-fi^^g)T/2] 
[2Tr(f-fjj^/2)T/2] 


(24) 


This  gives  the  range  bln,  or  location  of  the  target  In  rang'.  Its  center 
Is  slightly  displaced,  however,  due  to  the  Doppler  term  ^  f^.  In 
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Next  we  show 
the  s1nN»/s1na  factor 
of  (23)  in  figure  6 
for  a  moving  target. 
Notice  that  It 
repeats  itself  at 
the  pulse-repetition 
frequency ,  f^ ,  as 
one  would  expect  in 
any  pu1se-Dopp1er 
radar.  Note  that  this 
factor  contains 
only  Doppler  inform¬ 
ation  and  would  be 
identical  for  a 
target  at  any  range 
having  the  same 
radial  velocity. 

Finally  we 
show  the  product  of 
the  two  functions 
for  a  moving  target 
at  a  range  corre¬ 
sponding  to  fji^/2* 
Notice  that  the 
slow.y  varying  range  ^ 
function  essentially 
envelope  modulates 
the  Doppler  function, 
now  isolating  the 
target  in  both  range 
and  Doppler.  The  ( 
range  resolution  Is 
essentially  the  width 


Envelope  for 
'^moving  target 


sin  Nx 

”  Sirt  X 


(b) 


2V| 

■C"'c 


(d) 


r  V  - 


Envelope  for 
'stationary  target 


sin  Nx 
sin  X 


for  non-moving  forget 


fr  *1 

I  Spectrum  of 
(c)  repeated 


W  ''  v  ''f 

mf  -^iVRangebln  center 

•  for  non -moving  forget 

-►I  ^  »  Doppler  shift 

n[  11  for  target  in  m-th 
I  I  range  bin 


Figure  6.  Piatcrial  deeoription  of  long  single 
FFT  proaeeeing  over  N  pulaee. 


ifTT.' J  '<M;y#.-» j  ■■ '  i^y.- , 


. . . 


Of  the  broader  function,  corresponding  to  fj,  ■  1/T^  Hz  (1  .e. ,  a  1.5  km 
bln  here),  vthlle  the  Doppler  resolution  Is  essentially  the  width  of  the 
narrower  line,  corresponding  to  1/NT,.  Hz,  the  coherent  Integration  time. 

The  frequency  axis  after  the  long  FFT  can  thus  be  broken  up  Into 
M/2  coarse  range  bins  of  width  1/T,.  Hz;  within  each  range  bln,  finer 
frequency  divisions  then  correspond  to  the  Doppler  spectrum  of  the  target. 
In  particular,  there  are  N  Doppler  bins  per  range  bln,  corresponding  to 
a  Doppler  resolution  of  1/NT,.  Hz.  It  seems  proper  therefore  to  center 
each  range  bln  on  a  zero-Doppler  line.  The  centers  of  each  range  bln— 
as  shown  In  figure  6— are  thus  located  at  mf,,(mult1ples  of  the  PRF)  along 
the  frequency  axl  s ,  and  extend  +  f ,./2  away  from  thi  s  central ,  zero- 
Doppler  position.  Thus  we  can  take  the  plot  along  the  positive  frequency 
axis  and  divide  It  Into  M/2  pieces,  each  centered  at  mfp  where 
0  <  m  <  M/2.  Each  piece  then  represents  the  Doppler  spectrum  of  an  Indi¬ 
vidual  range  bln.  Or,  we  can  have  the  computer  do  the  "dividing"  for  us, 
displaying  each  range  bln  however  we  choose.  For  example,  range  bins 
could  be  lined  up  behind  each  other,  closely  spaced,  to  give  a  3-d1mens1on- 
al  range-Ooppler-IntensIty  display.  Note  also  that  each  range  bln— and 
the  resulting  Doppler  spectrum  thus  obtained— Is  similar  to  the  Fourier 
transform  over  a  given  column  In  the  preceding  section;  both  are  range 
bins  containing  a  Doppler  spectrum  vlth  the  same  resolution  and  width. 
THEREFORE  THE  TWO  PROCESSING  TECHNIQUES  YIELD  IDENTICAL  RESULTS. 

A  little  thought  will  also  show  that  this  technique  will  work  for 
distributed  targets.  For  example.  If  we  have  many  targets  over  several 
range  bins  but  at  the  same  velocity,  we  will  effectively  have  several 
K(f)  functions  In  (23),  but  centered  on  slightly  different  positions. 

The  slnNx/sIn®  functions  for  the  Doppler  will  be  Identical.  Thus  In  effect 
the  target  at  a  given  Doppler  will  appear  In  several  range  bins,  as  It 
should,  but  at  the  same  discrete  velocity  In  each. 

8.  NUMBER  OF  COMPUTER  OPERATIONS  REQUIRED 

The  possibility  exists  with  present  day  computers— especially 
"mlnlcr-iputers"  of  the  NOVA  and  HP  2110/2115  variety— that  the  range- 
Doppler  processing  described  above  can  be  done  In  real  time.  Such 
processing  for  HF  radars  has  In  fact  been  done  digitally  in  real  time  by 


-'■  \vi  r«.  ‘  '•■■-  r-T*.i': ■'•'iv'i.'Af.'srT^^iHipi- . ^ 


several  groups  for  both  discrete  targets  and  sea  scatter,  using  no  more 
than  a  single  HP  2115  minicomputer.  To  ascertain  whether  such  Is  possible 
for  a  given  application,  we  must  know  the  number  and  size  of  the  digital 
words  to  be  stored  and  processed  per  second. 

The  FFT  process  Is  known  to  require  Llog^L  operations  for  a  linear 
array  of  L  numbers.  Let  us  first  analyze  the  total  number  of  operations 
required  by  the  double  FFT.  We  first  do  an  FFT  on  a  pulse,  using  M 
samples;  this  requires  MIogjM  operations.  Next  we  begin  transforming 
over  each  of  the  M/2  columns;  each  now  contains  a  real  and  Imaginary 
word  for  a  total  of  M  words.  With  N  elements  In  a  column,  Nlog2N 
operations  are  required  for  the  FFT  on  each  column.  For  M  column  words, 
this  gives  MNlog2N  operations.  Thus  the  sum  of  operations  required  In 
the  first  and  second  sets  of  FFT  processing  Is 

MNlog2M  +  MNlog2N  »MN(log2M  +  logsN)  -  MNlog2MN  (25) 

operatl ons . 

The  number  of  operations  required  In  the  single  long  FFT  Is  simple 
to  calculate.  With  N  pulses  and  M  samples  per  pulse,  we  have  MN 
total  samples  per  transform.  This  therefore  requires  MNlog^MN  operations. 
THIS  IS  IDENTICALLY  THE  SAME  NUMBER  AS  FOR  THE  DOUBLE  FFT! 

Normally  the  FFT  requires  that  the  number  of  samples  to  be  transformed 
be  an  Integer  power  of  2.  For  the  double  FFT  process  therefore,  both 
M  and  N  must  be  powers  of  2  (e.g..  256  and  128,  32  and  64,  etc.;  Just 
so  M  and  N  Individually  are  greater  than  the  number  required  by  the 
sampling  rate  and  Doppler  resolution).  For  the  single,  long  FFT,  the 
product  MN  must  be  a  power  of  two,  and  hence  again  M  and  N  must 
Individually  be  powers  of  two. 

In  both  cases,  MN  elements  must  be  accumulated  and  stored  for 
processing;  this  dictates  the  size  of  the  r'SquIred  core  und/or  disc 
storage.  The  entire  number  of  MNlog  MN  operations  must  be  performed 
every  NT^  seconds  If  the  process  Is  to  be  done  1ii  real  time.  This 
requires  that  (Mlog2MN)/T^  computer  operations  per  second  be  done  (not 
Including  time  for  buffering  and  display  functions).  Thus  the  obvious 
way  to  reduce  the  require'd  data  rate— If  such  Is  necessary— Is  to  lower 
M,  the  number  of  range  bins.  Since  M  Is  equal  to  2f„B,vTy.,  we  must 
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reduce  f  >  the  maximum  IF  frequency  per  pulse.  This  does  not 

ll&flLX 

necessarily  require  one  to  reduce  the  range  resolution.  For  example, 

suppose  for  our  example  that  Instead  of  observing  a11  ranges  from  0  to 

150  km  with  a  1.5  km  resolution  (giving  M>200),  we  decided  that  we  only 

wanted  to  observe  the  window  between  126  km  and  150  km,  but  still  with 

1.5  km  resolution.  This  gives  conceptually  M»32  or  M/2-16  range  bins. 

To  achieve  this,  one  merely  slides  the  linear  sweep  delay  In  the  receiver 

so  that  Instead  of  varying  between  84  and  f  =100  Hz,  f  now  runs 

imax  1 

between  0  and  f  =16  Hz.  Then  the  M=32  samples  are  adequate  for 
max 

the  T^»l  second  pulse  repetition  Interval. 

Finally,  the  number  of  bits  required  per  word  also  affects  the  data 
rate  to  some  extent.  The  processor  dynamic  range  depends  upon  the  bits 
per  word  because  of  quantization  error.  Thus  the  dynamic  range  is  optimally 
6b  decibels,  where  b  is  the  number  of  (binary)  bits  per  word.  Currently 
about  80  dB  dynamic  range  can  be  realized  by  digital  processors  without 
too  much  difficulty,  requiring  14  bit  words  and  a  14-bit  A/D  convertor. 

9.  WINDOWING  AND  WEIGHTING 

In  all  of  the  preceding  sections,  we  assumed  a  square  pulse  at 
frequency  fj,  and  N  such  pulses  all  with  the  same  amplitude.  As  a  result 
we  arrived  at  sinx/x  and  sInNx/sinx  functions  in  the  frequency  domain 
for  the  target  echoes,  doth  functions  have  rather  high,  objectionable 
sidelobes:  the  first  sidelobe  of  the  sinx/x  function  Is  only  13  dB  down 
from  the  main  lobe,  while  the  average  sidelobe  level  of  the  sinNx/sinx 
function  between  main  lobes  Is  only  down  20  dB.  Thus  some  of  the  side- 
lobes  from  a  single  target— as  Illustrated  in  figure  6— are  quite  high 
and  could  be  mistaken  for  other  targets. 

The  remedy  for  this  is  the  same  as  that  taken  by  antenna  designers 
to  reduce  sidelobes:  use  an  amplitude  taper  across  the  original  function 
before  Fourier  transforming.  This  technique  is  currently  being  used  in 
nearly  all  radar  digital  processing  schemes.  The  common  amplitude 
taper— or  weighting— used  across  the  time  window  is  the  Taylor  weight 
(although  Hamming  and  cosine-squared  weights  (Blackman, 1958;  Nathanson, 

1969)  are  sometimes  used).  This  results  in  average  sidelobes  down  40-50  dB 
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below  the  main  lobe.  The  only  bad  effects  of  such  weight  ng  are  the  slight 
broadening  of  the  main  lobe  (by  as  much  as  401(  In  some  cases  at  the  3-d6 
point)  and  a  drop  of  1-2  dB  In  signal-to-nolse  ratio  due  to  attenuation  of 
the  original  received  signal  at  the  edges  of  the  wlndov^. 

For  both  types  of  processing  described  above,  two  weighting  functions 
are  normally  performed  digitally.  The  first  Is  to  weight  the  M  samples 
within  the  pulse  according  to  the  selected  function  (e.g.,  Taylor  weight¬ 
ing).  The  next  Is  to  weight  the  N  pulses  to  be  used  In  the  coherent 
Integration  by  the  selected  technique.  Both  weighting  processes  across 
the  two  respective  windows  of  T^  and  NT,,  seconds  are  normally  required 
to  keep  both  the  range  and  Doppler  sidelobes  unobjectionable. 

10.  RULES  FOR  SIGNAL  DESIGN 

Here  we  give  a  simple,  stepwise  procedure  for  calculating  the  signal 
parameters  required  for  a  given  set  of  backscatter  radar  or  sounder  specifi¬ 
cations.  We  assume  that  the  following  parameters  describing  the  system  are 
given:  (1)  f^..  the  carrier  frequency.  In  Hertz;  (11)  the  range 

window  width  to  be  calculated  and  displayed.  In  meters;  (111)  V|^,  the 
maximum  target  velocity  In  m/s;  (iv)  AR,  the  range  resolution  desired, 
in  meters;  (v)  Av,  the  velocity  resolution  desired.  In  m/s. 

With  these  parameters  given,  the  following  four  steps  are  to  be  used 
to  calculate  the  following  four  FM/CW  signal  and  processing  parameters: 

(1)  B,  the  signal  bandwidth,  or  frequency  excursion.  In  Hertz;  (11)  T, , 
the  pulse  repetition  Interval,  In  seconds;  (111)  N,  the  number  of  pulses 
of  period  T^  needed  for  a  single  coherent  processing  operation;  and 
(Iv)  M,  the  number  of  samples  needed  per  pulse  Interval,  T^. 

(1)  B  •  c/(2AR),  where  c  Is  the  wave  propagation  velocity  In  the 
medium. 

(2)  T^  ■  1/f,.,  where  f^  ■  fp^  being  the  maximum  target 

Doppler  shift,  given  by  fp,^,  ■  (2vn,/c)f^. 

(3)  N  »  Tg/T^,  where  T^.,  the  total  coherent  Integration  time  Is 
the  reciprocal  of  the  desired  Doppler  resolution,  Afp,  where 

Afp  ■  (2Av/c)fj.. 


(4)  M  •  2R^MR  samples  per  pulse  Interval,  T^,. 

In  the  above,  we  have  assumed  that  f^,  R^,  V|^,  AR,  and  Av  were 
a11  given  and  that  B,  T^,  N,  and  M  were  to  be  found.  In  practice,  the 
size  of  the  computer  and  data  handling  rate  will  often  limit  M  and  N. 
Thus  one  usually  Iterates  until  an  acceptable  compromise  Is  achieved,  i.e., 
he  varies  his  requirements  for  R^^,  AR,  and  Av  until  values  of  M  and 
N  are  obtained  within  the  capacity  of  his  machine. 

11.  SATISFACTION  OF  REQUIRED  ASSUMPTIONS 

In  the  course  of  the  analysis  herein,  certain  assumptions  were  made, 
upon  which  the  desired  output  Is  dependent.  If  these  are  not  satisfied, 
quadratic  and  other  types  of  distortions  will  result  which  reduce  or  limit 
the  achievable  signal -to-nolse  ratio.  Having  derived  B,  T^,  M,  and  N 
from  the  rules  of  the  preceding  sections,  one  can  quickly  check  the  follow¬ 
ing  criteria  to  see  whether  the  optimum  processing  gain  will  be  realized. 

(1)  }bT^(2v„/c)2n*  «  1, 

B(2v„/c)(2R^/c)N  «  1, 

B(2v„/c)T^/4  «  1 

Satisfaction  of  the  above  conditions  was  assumed  in  going  from  (12) 
to  (13)  for  the  phase;  If  one  or  more  of  these  conditions  are  not 
satisfied,  distortion  will  reduce  the  achievable  processing  gain. 

(2)  v^NT^  <  AR. 

This  merely  means  that  the  target  Is  not  traveling  so  fast  that  It 
moves  through  several  range  bins  within  one  coherent  Integration  period, 
NT^.  If  the  Inequality  falls.  It  simply  means  that  the  echo  will  appear 
In  several  range  bins,  but  with  a  proportionately  reduced  amplitude  in 
each. 

(3)  dv/dt  NT^  <  Av. 
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This  assumption— htratofort  unmantlonad— concorni  tht  rata  of 
change  of  radial  target  velocity  (or  radial  accelaratlon)*  It  has  been 
assumed  throughout  the  analysis  that  the  targets  under  consideration  have 
a  constant,  nonaccelerating  velocity.  Small  radial  accelerations  can  be 
tolerated,  but  If  dv/dt  Is  sufficiently  large  that  the  above  Inequality 
falls,  then  the  echo  will  appear  spread  Into  several  Doppler  bins  with 
proportionately  reduced  amplitude  In  each. 

12.  SUMMARY 

Despite  statements  often  seen  concerning  "chirp"  (I.e.,  linearly 
swept  frequency)  signals  used  with  microwave  radars,  there  te  no  anbiguity 
between  target  range  and  velocity  for  processing  done  In  the  straight¬ 
forward  digital  manner  described  In  this  note.*  Furthermore,  two  seemingly 
different  digital  processing  schemes  are  described  and  analyzed  herein, 
which  will  produce  exactly  the  same  pulse-Ooppler  (range-velocity)  output. 
Both  employ  the  FFT;  the  first  uses  a  shorter  FFT  many  times,  while  the 
second  uses  only  one  long  FFT  to  produce  the  same  coherent  pulse-Doppler 
map.  Both  techniques  work  equally  well  for  discrete  targets  (such  as  an 
aircraft),  as  well  as  for  continuous  or  distributed  target  complexes 
(such  as  ocean  waves,  rain,  atmospheric  tu^'bulence,  etc.),  and  display  the 
targets  In  their  appropriate  range-velocity  perspective. 

Identically  the  same  total  number  of  FFT  operations  is  required  for 
both  techniques;  the  same  data  rate  (A/D  convertor  rate)  Is  required  In 
each  case  also,  1.e.,  (2R^/AR)x2x(2vn|/c)fj.  words  per  second.  Here,  R^ 

Is  the  range  window  length  to  be  examined,  AR  is  the  range  resolution 
desired,  v^^  Is  the  maximum  target  velocity  to  be  encountered,  c  Is  the 
free  space  wave  propagation  velocity,  and  f^.  Is  the  carrier  frequency. 

The  choice  of  whether  to  use  the  multiple  vs  the  single  FFT  processing 
technique  then  rests  with  the  availability  of  appropriate  equipment.  For 

*Perhap8  the  ambiguity  occurring  In  the  microwave  aystema  is  attributable 
to  the  analog  pulse  compression  techniques  commonly  employed  there,  such 
as  the  dispersive  delay  line„  Here  the  technique  used  is  more  properly 
described  as  a  coherent  correlator  followed  by  pulse-Doppler  processing, 
rather  than  time-domain  pulse  compression.  The  difference  between  the 
two  techniques  results  in  the  elimination  of  the  ambiguity  for  type  of 
processing  described  here. 
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•xtinplti  tiMl)  cwpMttri  iHy  Ut  HmUtd  In  tht  ilit  of  «  ttinglo  FfT  thty 
cm  hmdlti  In  this  ciiOi  the  muUlplo  fFT  itchnlqut  htvino  imilUr 
unit  tlio  wey  bo  roqulrtd.  On  the  other  hind,  iptclil  hird*w1rtd  FH 
coMputort  aro  currently  •vellable  (called  ^FFT  boxea'*),  Theie  cm  perform 
a  fairly  largoi  fixed-length  traniforix  very  rapidly  bacauie  of  their 
ipeclaMied  conitructloni  and  are  uaed  at  one  cofapnnent  In  the  overall 
digital  procetting  lyitem.  Here,  the  tingle  long  FF1  It  utually  More 
efflelent  beautt  the  need  for  continual,  Interactive  itorage/retrleval 
of  eloMenti  In  matrlx/faihlon  demanded  by  the  multiple  FFT  tcheme  It 
eliminated. 
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